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1. Introduction 

Precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements provides 
a stringent test of the Standard Model to search for new physics. One of the elements |V,, 4 | can be 
determined from the K^nlv decays through a theoretical calculation of the normalization of the 
vector form factor /+(0), which parametrizes the K — > % matrix element 

{nip'W^Kip)) = (p + p^^U^ + ip-p')^^) (q 2 = (p-p') 2 ). (1.1) 

In this article, we report on our calculation of / + (0) in Nf = 2+ 1 lattice QCD within 1 % accuracy. 
In order to demonstrate the reliability of this precise calculation, we also examine the consistency of 
other quantities, such as /_ (0) and the form factors' shape, with chiral perturbation theory (ChPT) 
and experiments. 



2. Calculation of form factors 

We employ the overlap quark action to exactly preserve chiral symmetry for a straightforward 
comparison of f{ + _}(q 2 ) with ChPT. Numerical simulations are accelerated by modifying the 
Iwasaki gauge action [1] and by simulating the trivial topological sector [1, 2]. Effects of the 
fixed topology are suppressed by the inverse of the lattice volume x N t [2] . We use our gauge 
ensembles generated at a single lattice spacing a = 0. 112(1) fm, which is determined from the Q. 
baryon mass, and at a strange quark mass m s = 0.080, which is close to its physical value m s . p h ys = 
0.081. We simulate four values of degenerate up and down quark masses m uc i = 0.015, 0.025, 
0.035 and 0.050 that cover a range of the pion mass 290-540 MeV. At each m lu i, a lattice size of 
16 3 x48 or 24 3 x48 is chosen to control finite volume effects by satisfying a condition M K L>4. 
The statistics are 2,500 HMC trajectories at each combination of m m i and m s . 

We calculate the scalar form factor fo(q 2 ) =f+(q 2 ) + f-{q 2 )q 2 / {M 2 * — Ml) at q 2 = q 2 asx = 
(M K -M„) 2 , f + (q 2 ) and ${q 2 ) =f-(q 2 )/f+(q 2 ) at q 2 < q 2 max from the following ratios [3] 

Cf*(A^';0,0)C; g (Af,Af';0,0) {M K + M n ) 2 2 2 

Cf K (At,At';0,0)C^(At,At';0,0) a,,at'^ 4M K M n /oWmaxJ ' {ZA) 

C£*(Af,A^p,p')C*(Af,0)C*(A?',0) J E K +E' n E K -E'. 2 .\ f + {q 2 ) 



Cf*(At,At';0,0)C K (At,p)C*(At',p') \M K +M K M K + M n hKq ' J f {q 2 V 



imaxy 



= C^(At,At';p,p')Cf K (At,At';p,p') ) E K + E' K { p + p ') k + {p- p ') k ^ q 2 ) 
k Cf*(At,At';p,p')C§ K (At,At';p,p') (p + p') k E K + E' n + (E K - E'^q 2 ) 

where E^} (P = n or K) represents the meson energy with the spatial momentum pW . Note that 
these observables are sufficient to determine // + _m. (q 2 ) at simulated values of q 2 , except at <7m ax , 
where R and R^ have no sensitivity to %(q 2 ). Two- and three-point functions are defined as 

C p (At,p) = J_£^(^(x',r + A f )^(x,0)^' p(x '- x) ! (2.4) 



C PQ (At,At';p,p') = -i-££(^ e (x /, ,? + A? + AO^(x / ,f + A0^(x,0> 

xe - ( p'(x"-x')-<p(x'-x) ) (2 _ 5) 
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M ~ 540MeV, 9=1.60, 6' = 0.00 



l_ M K ~ 540 MeV, m ( ~ m ? phys , mj - m s ~ - 25 MeV 
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Figure 1: Left panel: Monte Carlo history of normalized reweighting factor w(m' s ,m s ) at our largest m U( j. 
Four lines show data with different numbers of noise samples (A^,-). Right panel: statistical fluctuation of 
Cf K and R before (open symbols) and after (filled symbols) reweighting. W represents the twist angle used 
to induce the initial (final) meson momentum p". 



where ffp^(x,t) (P,Q = n orK) represents the meson interpolating operator £ r 0(r)g(x+r,?)y5g'(x,?) 
with an exponential smearing function 0(r) =e~°- 4 l r L We use the all-to-all propagator [4, 5] to re- 
markably improve the statistical accuracy of the correlators and, hence, form factors [6, 7]. 

In order to precisely determine /+ (0) , we explore the most important kinematical region q 2 ~ 
by using the twisted boundary conditions (TBCs) [8] 

q(x+N s k,t)=e w q(x,t), q(x + N s k,t) = e~ w q(x,t) (k = 1,2,3), (2.6) 

where k is a unit vector in the k-th direction. For simplicity, we use a common twist angle 6 in 
all three spatial directions and take four (three) values of 8 (including 6=0) to cover —0.1 < 
(? 2 [GeV 2 ] <^ ax on the 16 3 x48 (24 3 x48) lattice. 



3. Reweighting 

The strange quark mass dependence of the form factors is studied by simulating a different 
value m' s = 0.060, which is about 25 MeV smaller than m s = 0.080. We employ the reweighting 
technique in which an observable & at m' s is calculated on the gauge configurations at m s as 

, ^ i ^ / / nv / / \ w(m'm s ) . , , det[D(m')l ,„ .. 

(0) m , = {0w(m! s ,m s ms , w(m' s ,m s ) = \" , w m>, = ' , (3.1) 

(w(m' s ,m s )) ms det[D(m,)J 

where (• • • )^(,) and D(m.['') represent the Monte Carlo average and the overlap-Dirac operator at 

the strange quark mass mi'\ respectively. We consider decomposing the reweighting factor w 
into contributions of low- and high-modes of D as w = wi ow Whigh- The low-mode contribution 
wiow = Ilk ^k{m' s ) / Hk ;^k( m s) ls exactly calculated by using 160 (240) low-lying eigenvalues Xk 
on 16 3 x48 (24 3 x48). A noisy estimator is employed for wg igh = (l/A^lf^expf-^c;,.) 1 "^- 

\){PE, r )/2] with Q. = D(m s )^D(m' s )- i y Dim',,)- 1 D(m s ), where {<j;i ,...,£#,} is a set of Gaussian 
noise vectors, and P projects them to the eigenspace spanned by the high-modes. 
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The left panel of Fig. 1 shows the Monte Carlo history of the normalized reweighting factor 
VP, which appears in the expression of the observable in Eq. (3.1). With our simulation setup, VP has 
small dependence on the number of the noise vectors N r . From this observation, we set N r = 10 to 
precisely calculate VP on each gauge configuration. 

The right panel of Fig. 1 shows how reweighting affects the statistical accuracy of the correlator 
Cf 71 and the ratio R. The fluctuation of Cf 71 is largely enhanced by VP, which is typically in the 
range of [0.5,2.0]. We observe, however, that the enhanced fluctuations are largely canceled in 
the ratio R. Consequently, the statistical accuracies of the form factors are not largely impaired 
by reweighting: typically < 1.0 % (20 %) for f{ + fi}(q ) (<§ (q 2 )) before reweighting, and < 1.5 % 
(30 %) after reweighting. 



4. q 2 dependence of form factors 

Our results for fo{q 2 ) at the smallest m U( j are plotted as a function of q 2 in the left panel of 
Fig. 2. We use TBCs to simulate small values of \q 2 \, where contributions of higher orders in q 2 
are small and our data are well described by any of the following parametrization forms 

/o(^ 2 )=/o(0)(l+c^ 2 ), Mq 2 )=M0)(l+ Cl q 2 +c 2 q 4 ), f (q 2 ) = f °}^ . (4.1) 

A different form /+ (0) { 1 / ( 1 - q 2 /M^ ) + ci<7 2 } with the K* pole plus a polynomial correction is 
also tested for f+{q 2 ). We confirm a good agreement among /+(0) and /o(0) obtained from these 
parametrizations as shown in the right panel of Fig. 2. In this report, we employ a simultaneous fit 
using the quadratic form for fo and the form with the K* pole for /+ to determine the normalization 
/+(0) and its slope /+(0). The uncertainty of this interpolation is estimated by using the different 
parametrization forms and turns out to be similar to or smaller than the statistical error. 




Figure 2: Left panel: scalar form factor fo{q 2 ) at our smallest m uc / as a function of q 2 . We plot interpola- 
tions listed in Eq. (4.1) together with /o(0) determined from the quadratic parametrization (diamond). Our 
result on a smaller lattice 16 3 x48 is also plotted to examine finite volume effects (triangle). Right panel: 
comparison of /+(0) (= /o(0)) obtained by using different parametrization forms. 
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In the left panel of Fig. 2, we also plot a result of fo(q 2 ) obtained on a smaller lattice 16 3 x48. 
The difference from 24 3 x48 can be attributed to the conventional finite volume effect as well as 
the fixed topology effect, but turns out to be insignificant (0.8 %, 2.1 a). The finite volume effect 
remaining on the larger volume 24 3 x48 is estimated as < 0.3 % by assuming the 1 /N^N t scaling 
of the fixed topology effect and can be safely neglected in the following analysis. 



5. Chiral extrapolation of /+ (0) 

There are two popular choices of the expansion parameter in ChPT: t,p = Mpj{A%F % ) 2 and 
Mp/(47iFo) 2 (P = n,K,ri), where F$ is the decay constant in the St/ (3) chiral limit. The latter 
choice largely enhances the chiral corrections in / + (0) [7] and leads to worse convergences in the 
expansions of Mp and fp [9]. We therefore use the former and denote the chiral expansion as 
/+(0) = 1 +/2 + Af, where fi and Af are 0(t,p) and higher order contributions, respectively. 

The Ademollo-Gatto theorem /+(0) — 1 °< (m s — m u d) 2 [10] guarantees that fi [1 1] 



- ^ ( H K7t + HfCri ) , Hpq — - ^ P g ^ Q [ 1 + p_2 P ^2 ln 



op bQ 



(5.1) 



is written in terms of physical observables ^i^k,^} an d does not contain low-energy constants 
(LECs) in the ChPT Lagrangian. The chiral expansion of /+(0) is then nothing but the parametriza- 
tion of Af. In the left panel of Fig. 3, we examine the quark mass dependence of Af divided by 
[M\ — M 2 ) 2 which is motivated by the Ademollo-Gatto theorem. The mild dependence suggests 
that our data can be described by a simple constant fit Af j {M^—M 2 ) 2 = C4 which gives rise to the 
0{£,p) analytic term in /+(0). We also confirm that the chiral extrapolation of A//(Mf — M 2 ) 2 is 
not largely modified by including the following 0{E,p) logarithmic and 0(%p) analytic corrections 

A/ '/ '(Ml - Ml) 2 - c 4 = c' A n \og[£, K ], or c^log 2 ^], or c 6 ^ K , or c^^+c^k- (5.2) 

In this report, we employ a parametrization Af j{M\ — M 2 ) = C4 + c^,^^ K , all parameters of 
which are determined reasonably well with an acceptable value of # 2 /d.o.f. ~ 1.6. The systematic 




Figure 3: Chiral extrapolations of A// (M\ —M 2 ,) 2 (left panel) and /+(0) (right panel). Circles and squares 
show data at different values of m s . In the right panel, we also plot /+(0) from recent calculations in 
N f =2+1 [12] andjV/ = 2 [13] QCD. 
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uncertainty of this extrapolation is estimated as the largest deviation in / + (0) among the different 
parametrization forms discussed above. The discretization error in the SU (3) breaking effect f2+ 
Af is estimated by an order counting 0((<2Aqcd) 2 ) with Aqcd ~ 500. We then obtain 

/+(0) = 0.959(6) stat (4) chiral (3)^ , \V„\ = 0.2256(19) theory+eX p't, (5.3) 

where we use |V Mi / + (0)| =0.2163(5) determined from the K^nlv decay rates [14]. Note that 
previous calculations in Nj ■ = 2 + 1 [12] and Nf=2 [13] QCD are consistent with our result. 

6. Comparison of {r 2 )y n and E, (0) with ChPT and experiments 

In order to demonstrate the reliability of the 1 % calculation of /+ (0) , we compare our data 
of the normalized slope (r 2 )f K = 6/^(0) / f + (0) and £(0) with ChPT and experiments. In contrast 
to /+(0), the Ademollo-Gatto theorem is not applicable to these quantities, and unknown LECs 
appear already in their leading chiral corrections. For (r 2 )y K , we parametrize its higher order 
contributions by simple analytic terms to avoid unstable chiral extrapolations 

(r*)f = l2L r 9 /Fl + "chiral logarithms" + d„^ + d K &. (6. 1) 

We refer to Ref. [11] for the explicit expression of the chiral logarithmic terms. The left panel 
of Fig. 4 shows that this form describes our data reasonably well. It also suggests that the chiral 
behavior of (r 2 )y n is significantly modified in our simulation region 290<M^[MeV] < 540 by 
the 0{t,p) contribution at two-loop order in ChPT. We note that significant two-loop corrections 
have been also observed in our studies of the pion and kaon charge radii [5, 6]. We confirm a 
good agreement of the extrapolated value of {r 2 )y n with experiment [15]. Our fit result L' 9 x 10 3 = 
4.3(0.6) s tat(0.3) sys is also consistent with the phenomenological estimate 5.9(0.4) [16]. 

In this report, we parametrize the quark mass dependence of B, (0) by a simple linear form 
% (0) = do + d\ (Mg — M 2 ), which is motivated from the ChPT expression of the leading analytic 
terms oc m\ —M 2 [11]. This form describes our data reasonably well as seen in the right panel of 
Fig. 4. We obtain do = —0.022(25) confirming that <^ (0) vanishes in the SU (3) symmetric limit, as 




[GeV 2 ] M K~ M n [° ev2 ] 



Figure 4: Chiral extrapolation of (r 2 )y K (left panel) and |(0) (right panel). In both panels, the diamonds 
represent our result extrapolated to the physical quark masses m uc [ and m s , which should be compared to the 
experimental values plotted by the stars. 
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expected. The extrapolation to the physical point yields £ (0) = —0.095(5) which is consistent with 
the experimental value —0.125(23) [17]. 

7. Summary 

In this article, we report on our calculation of the kaon semileptonic form factors. The nor- 
malization /+(0) is calculated within 1 % accuracy by utilizing the all-to-all quark propagator, 
reweighting and TBCs. The reliability of this precise calculation is checked by confirming a good 
consistency of {r 2 )y n and <§(0) with experimental results. 

Since we observe significant two-loop contributions in the chiral expansions of /+(0) and 
{r 2 )y n , it is interesting to apply two-loop ChPT formulae to our data. To this end, our use of 
the overlap quark action is advantageous, since exact chiral symmetry enables us to use the two- 
loop formulae in the continuum limit without any additional terms at finite lattice spacings. This 
provides a theoretically clean comparison between lattice QCD and ChPT at the higher order. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 1 1-05) as well as on Hitachi SR16000 at YITP in Kyoto University. This 
work is supported in part by the Grants-in-Aid for Scientific Research (No. 21674002, 21684013), 
the Grant-in- Aid for Scientific Research on Innovative Areas (No. 2004: 20105001, 20105002, 
20105003, 20105005, 23105710), and SPIRE (Strategic Program for Innovative Research). 
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